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We study classical hard-core dimer models on the square lattice with links extending beyond 
nearest-neighbors. Numerically, using a directed-loop Monte Carlo algorithm, we find that, in the 
presence of longer dimers preserving the bipartite graph structure, algebraic correlations persist. 
While the confinement exponent for monomers drifts, the leading decay of dimer correlations remains 
1/r 2 , although the logarithmic peaks present in the dimer structure factor of the nearest-neighbour 
model vanish. By contrast, an arbitrarily small fraction of next-nearest-neighbor dimers leads to 
the onset of exponential dimer correlations and deconfinement. We discuss these results in the 
framework of effective theories, and provide an approximate but accurate analytical expression for 
the dimer correlations. 
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I. INTRODUCTION 

Dimer models have a venerable history in statistical 
physicsi2i2ii More recently, they have also emerged as 
central models in modern theories of strongly correlated 
quantum systems, e.g., high-temperature superconduc- 
tors and frustrated antiferromagnets^SiZ&SiiSiiiiiSii&M 
There, dimers can represent singlet-forming electron 
pairs. In order to model the quantum fluctuations of 
the dimers and realize a short-range version of Ander- 
son's resonating valence bond (RVB) stateji£*i& Kivelson, 
Rokhsar, and Sethna introduced a Hamiltonian with a 
resonance term which flips pairs of parallel dimers on the 
two-dimensional (2D) square latticed The purely clas- 
sical dimer model also retains it relevance here; it was 
shown that the equal-weight sum over all dimer configu- 
rations is a ground state of the Hamiltonian when the res- 
onance strength equals the potential energy cost of each 
resonating pair of dimers (the RK point) & The equal- 
time dimer correlations of the quantum dimer model are 
thus those of the classical dimers jii and the imaginary 
time dynamics can be related to a classical random walk 
in configuration space itLlSLSiL 

The dimer pair correlations of the classical square lat- 
tice model decay with distance as r~ 2 and two inserted 
test monomers are correlated with each other as r -1 / 2 . 4 
As it turned out, away from the RK point the dimers form 
long-range order and the monomers are exponentially 
confinediiSi^ Hence this system does not give rise to 
the desired RVB state with no broken lattice symmetries 
and deconfined monomers (corresponding to spin-charge 
separation^). More recently, it was shown that a true 
extended RVB phase does appear in the quantum dimer 
model on the triangular latticed Following this insight, 
further work has been carried out in order to characterize 
classical and quantum dimer models on various lattices. 
Moreover, there are currently intense activities in gauge 
theories related to quantum dimer models&iSii^ii 2 ^ 



To date, research on dimer models in two dimensions 
has focused mainly on planar lattices, i.e., ones that have 
no intersecting links. This class of models can be analyti- 
cally solved, using Fermionic path integrals, with the aid 
of a theorem by Kastelyn^ Thus the quantum ground 
states at the corresponding RK points are characterized 
as well. This body of analytical work has established a 
number of properties of classical dimer models which are 
believed to hold true even for non-planar lattices. The 
basis for this belief is provided by effective theories in- 
corporating those exact results. Such effective theories 
are expected to be more robust than the exact solubility, 
and therefore not crucially dependent on lattice planarity. 
Consequently, there have been relatively few numerical 
studies testing these beliefs^ although it is not neces- 
sarily the case that the range of behaviours unearthed 
so far in fact exhausts all possibilities. In this paper, we 
provide a detailed study addressing some of these issues. 
We do this by considering two extensions of the square 
lattice dimer model which preserve its full symmetry. 

The first case we consider here is a nonbipartite model 
with nearest-neighbour (Ni) and next-nearest- neighbour 
(N2) dimers. Such bonds will inevitably appear in re- 
alistic systems away from the limiting cases represented 
by the nearest-neighbor (N±) dimer models. Introducing 
next-nearest-neighbor links along only one of the diag- 
onals destroys the square lattice symmetry and instead 
yields the topology of the triangular one. It has already 
been shown that an arbitrarily small fraction of such di- 
agonal dimers destroys the critical Ni square-lattice state 
and leads to deconfinement ^2^^ as in the isotropic tri- 
angular latticed The model with links along both diago- 
nal directions, i.e., the full 2D square lattice with N\ and 
N2 bonds, is not solvable by Kastelyn's theorem. 3 One 
would suspect (see e.g. Ref. |28j) that the critical state is 
immediately destroyed in this case as well but this has 
not been explicitly demonstrated. Similarly, one might 
speculate that introducing longer bonds between the two 



sublattices does not lead to deconfinement, but there are 
no exact results to back this up. For this reason, we also 
study the bipartite lattice with Ni and fourth-nearest- 
neighbors (N4, of which there are eight per site). 

In both the above cases cases, we use an efficient 
directed-loopSiSi Monte Carlo algorithm to sample the 
full space of hard-core dimer configurations, with fugaci- 
ties Wi assigned to the different dimer types Ni . We then 
discuss our results in the framework of effective and sim- 
ple microscopic analytic treatments. We confirm that 
a low (most likely infinitesimal) 3 " concentration of N2 
dimers leads to deconfinement; in particular, we find that 
the finite correlation lengths for monomers and dimers in- 
duced by the presence of the N 2 dimers scales inversely 
to their fugacity; £ oc \jwi- By contrast the presence of 
N4 dimers preserves the algebraic decay of dimer correla- 
tions, ^1/r 2 . These correlations have a complex dipolar 
structure in reciprocal space, which implies that they do 
not lead to a logarithmic divergence of the structure fac- 
tor for W4 7^ 0; such a divergence is present for the pure 
N\ model. For these observations, we provide an approx- 
imate analytical theory based on the height model in the 
Coloumb gas representation and a self-consistent mean 
field solution of a large-TV generalization of the dimer 
model. 

Our numerical simulations show that the confinement 
exponent governing the monomer correlation function 
drifts, changing continuously from —1/2 in the pure N\ 
case to —1/9 (to high numerical accuracy) for the pure 
A4 model. In the language of height models, this cor- 
responds to a decrease of the stiffness induced by the 
addition of longer-ranged dimers. Note that for the pure 
JV4 model, like for the soluble dimer and six-vertex mod- 
els, the stiffness appears to be a simple rational number. 
The efficient directed loop algorithm provides an efficient 
direct determination of this quantity in general. 

The outline of the remainder of the paper is the fol- 
lowing: In Sec. II we introduce the directed-loop al- 
gorithm for a wide class of dimer models, and in in 
Sec. Ill we present simulation results. In Sec. IV we 
qualitatively explain the numerical findings based on 
the height-model/Coulomb-gas representation of bipar- 
tite dimer models and their large-N generalizations. We 
conclude in Sec. V with a summary and a brief discussion 
of the relevance of our study to RVB physics. 



II. DIRECTED-LOOP ALGORITHM 

This algorithm is an adaptation of a quantum Monte 
Carlo methodf^i 3 ^ with the same name, in which updates 
of the system degrees of freedom are carried out along a 
self-intersecting path at the endpoints of which there are 
defects not allowed in the configuration space contribut- 
ing to the partition function. When the two defects meet 
they annihilate, the loop closes, and a new allowed config- 
uration is created. The conditions for detailed balance in 
the process of stochastically moving one of the defects are 
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FIG. 1: (a) Labeling of the links of the N1-N2 dimer model, 
(b) A step in the directed-loop update, in which the entrance 
to the vertex is at site 3 and the exit is at site 8. The thick 
bond indicates the location of the dimer. 



expressed as a coupled set of directed-loop equations, 29 
The applicability of this scheme to dimer models was 
first realized by Adams and Chandrasekharani^I Here an 
algorithm for the multi-length dimer problem will be pre- 
sented; a simplifying representation of the dimer config- 
urations will also be introduced. The algorithm will be 
described only for the case of the N1-N2 model, but the 
scheme applies directly to any range of the links. 

The links connected to a given site are labeled as shown 
in Fig. ^a). For a dimer configuration on a periodic 
L x L lattice, each site is numbered according to the type 
of dimer it is connected to. For any given site, the other 
member of the same dimer can then be easily found. The 
central object in the directed-loop algorithm is a vertex, 
which in this case consists of a site and all its surrounding 
sites to which it can be coupled by a dimer. The state 
of the vertex is the number (1-8) assigned to the central 
site. A step in the directed-loop algorithm is illustrated 
in Fig. ^b). The vertex is entered through the dimer, 
and one of the seven other surrounding sites is chosen as 
the exit. The dimer is then flipped from the entrance to 
the exit, the central vertex site remaining connected to 
it. The exit site already has another dimer connected to 
it, and this site now becomes the entrance in the next 
step of the algorithm. In each step a defect is hence 
moved one link ahead, leaving behind a healed dimer 
state. The first entrance site is chosen at random and the 
corresponding defect remains stationary. It is annihilated 
when the moving defect reaches it, whence a new allowed 
dimer configuration has been generated. In the present 
case the defects are monomers, and the intermediate two- 
monomer configurations can thus be used to determine 
the correlations between two inserted test monomers. 

The key to an efficient algorithm of this kind is that 
the probabilities for random selection of the seven possi- 
ble exit sites can be chosen in such a way that detailed 
balance is satisfied without any further accept/reject cri- 
terion. Each step then moves one dimer, and a full loop 
can accomplish very significant changes to the dimer con- 
figuration. The directed-loop equations 2 - give the con- 
ditions for detailed balance in terms of weights a.jk for 
the processes in which a vertex in state j is entered at 
site j and exited at k (transforming the vertex into state 
k). The actual probabilities Pjk = cijk/wj, which implies 



Detailed balance is satisfied if a 



jk ®kj • 



In the N1-N2 model, the vertices can be classified as even 
(e) or odd (o) according to the numbering of Fig. ^a); 
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there are then four weights: a ee , a 00 , a eo , a oe . In princi- 
ple, one can include "bounce" processes where j = k, but 
in the present case they can be excluded. Including only 
the seven no-bounce exits, the directed-loop equations 
reduce to 
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This system is underdetermined and has an infinite num- 
ber of postive-definite solutions. Here the following solu- 
tion will be used: For w\ > u>2, 
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while for W2 > W\, 



a Q e = Cleo = W2/7, 

( Wl - 4w 2 /7)/3, 
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(w 2 - 4«;i/7)/3. 
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There is no guarantee that this is the best solution, but 
the resulting algorithm performs very well and allowed 
for studies of lattices with ~T0 6 dimers. The fugacity of 
the A^-bonds is set to unity (except in the case of the 
pure N2- and A4-models) . The program was tested using 
known results for the pure ATi case^ and by comparing 
with local Metropolis simulations for small lattices. 
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FIG. 2: Dimer (a) and monomer (b) correlations along the 
direction (r, 0) for the N1-N2 model at W2 = e -4 , compared 
with those of the pure iVi model. The solid lines show the 
agreement with the known power-law- for the Ni model. Note 
that for the pure Ni model M(r) = for even r. The results 
were obtained using L — 1024 lattices. 



III. SIMULATION RESULTS 



Dimcr-dimer correlations in the full close-packed sys- 
tem and monomer-monomer correlations in the system 
with two test monomers will be discussed. The monomer 
correlations M(r) were obtained by accumulating the dis- 
tances between the stationary and the moving monomer 
in the directed-loop update. As has become customary, 
the normalization M(r = 1) = 1. Several types of dimcr- 
dimer correlations can be defined. Here -Ds(r) will be 
defined in the following way: If site i is connected to a 
dimer in the set E, a variable s(i) — 1, otherwise s(i) = 0. 
The correlation function Z?£ (ry ) is then given by the con- 
nected correlator (s(i)s(j)). Results will be presented for 
cases where E contains a single dimer or half of the dimers 
of a given type (i.e., two neighboring Ni or N2 dimers or 
four neighboring N4 dimers connected to a site i) For the 
N1-N2 model, the correlation functions D\, D2, Da, and 
Db are thus defined corresponding to the sets {1}, {2}, 
A = {1,3}, and B = {2,4}. Analogous definitions are 
used for correlations £>i, Da of N\ dimers and D4, Dp. 
of A4 dimers in the N1-N4 model. 

A. Deconfinement in the N1-N2 model 

In Fig. |2Ia), dimer correlations Da for the nonbipar- 
tite N1-N2 model with W2 — e~ 4 (corresponding to a con- 
centration P2 ~ 0.3% of N2 dimers) are compared with 



those of the pure Ni model. There is a clear deviation 
from the r~ 2 decay, showing that the behaviour of the 
N1-N2 model is fundamentally different even at this very 
low concentration of N2 dimers. As shown in Fig.[2Ib), 
the test monomer correlation approaches a nonzero con- 
stant, i.e., the system is deconfmed in contrast to the 
algebraically confined Ni model. The very significant 
changes seen already at a very low concentration of N2 
dimers suggest that an arbitrarily small concentration 
indeed causes deconfinement. 

How does deconfinement set in upon inclusion of N2 
dimers? To address this question, we have simulated the 
N1-N2 model for different (small) values of W2 ranging 
from e~ 3 w 1.5% to e~ 6 « 0.0086%. We analyze the 
correlation functions using a data collapse with a scaling 
function of the form r -Q $(r/£). With £ ~ 1/^2 and a 
the respective exponent for the correlations at W2 = 0, 
we obtain a reasonable data collapse, as shown in Fig. 
The data for monomers is easier to handle as the alge- 
braic decay is less rapid. We have also tried to improve 
the data collapse by adjusting the exponent a. We find 
the best-fit values a = 2.0 ± 0.1 and 0.55 ± 0.05 for the 
dimers and monomers, respectively, consistent with them 
remaining at their RK point values in this regime. 
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FIG. 3: Scaling collapse of monomer (top panels) and dimer 
(bottom panels) correlations as the fugacity W2 is varied. The 
two distinct sets of points in the upper panel correspond to 
even and odd distances. 



B. Correlations in the N1-N4 model 

Turning now to the bipartite N1-N4 model, its dimer 
correlations Da{t) (involving N\ dimers) and -Drj(r) (in- 
volving N4 dimers) are compared with Da(t) of the pure 
Ni model in Fig. 21 Both correlation functions appear 
asymptotically to decay according to the same form r~ 2 
as Da(t) of the pure N\ model. Interestingly, although 
Dr>(r) is much smaller in magnitude than Da{t) at small 
W4, it reaches the asymptotic form at shorter distances. 
In Fig. [3] the Fourier transforms of the correlation func- 
tions D\ and D4 (involving a single dimer) are shown 
for the pure N\ and N4 models. For the N\ model, the 
dominant correlations are at q = (tt, 0). Using finite-size 
scaling, a logarithmic divergence of the peak height with 
the system size can easily be observed (not shown here). 
As this feature masks some non-diverging correlations 
on the contour plot, we have cut it off to display more 
clearly, most notably, a non-analytic bow-tie feature at 

(7T,7r). 

In the N4 model, the logarithmic peak is absent, and 
the correlations display a rather more complex structure, 
exhibiting very broad peaks, and again a bow-tie visible 
at the zone corner; this is rotated with respect to the 
bow-tie for Ni dimers. These features show almost no 




FIG. 4: Various dimer correlations along the direction (r, 0) 
in the N1-N4 model (on L = 1024 lattices). The fugacities 
lt>4 = e -4 and e -3 correspond to concentrations P4 « 0.024 
and 0.065, respectively. The solid line shows the asymptotic 
form r~ 2 . 




FIG. 5: Fourier transform of the dimer-dimer correlation D\ 
of the Ni model [left, for dimers pointing along (1,0)] and 
D4 of the N4 model [right, for dimers pointing along (2, 1)] 
calculated on L = 128 lattices. In the left graph, the log- 
divergent peaks at (n, 0) and (it, 2tt) have been cut-off at the 
value 2.0 (for this lattice size the peak value is m 3.12) 



size dependence up to the largest size (L — 2048) studied, 
although one might naively have expected a logarithmic 
divergence here as well, considering the r~ 2 real-space 
correlations. 

Whereas the dimer correlations decay with a leading 
r~ 2 dependence, a drift in the exponent of the algebraic 
decay is manifest in the monomer correlations. Results 
for M(r) at r = L/2 — 1 are shown multiplied by L a in 
Fig.UJa). Here a is adjusted to give a flat L dependence, 
and hence the critical form M(r) ~ r~ Q is extracted. 
The known a = 1/2 is used for the pure N\ model. For 
the pure .ZV4 model the exponent is consistent with a = 
1/9 (to an accuracy of 1%). To show that the exponent 
changes from 1/2 already at a low concentration of 7V4 
bonds, results for W4 = e -5 (p 4 w 0.9%) are scaled with 
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FIG. 6: Finite-size scaling of monomer correlations in the N\- 
N4 model, (a) The best scaling for three different cases, (b) 
Failure of a scaling with the pure Ni exponent a — 1/2 for 
the N1-N4, model at low fugacity wa- 



a = 1/2 in Fig. Eft)). This scaling clearly fails, and is 
instead consistent with a w 0.485. 



IV. ANALYTICAL RESULTS 

We now present an analytical treatment of the numer- 
ical findings. We will in turn discuss exponents, cor- 
relation lengths, amplitudes and the angular form of 
the dimer correlations. Our discussion is based on the 
height model/Coulomb gas representation of bipartite 
dimer models on one handf^i^iS and on their large- N 
generalisation on the other m& The former is based on a 
mapping of the dimer configurations onto a fluctuating 
height surface and accounts for the gross features of the 
correlations described above. The latter reproduces in 
detail the features of the correlations of the N4 model. 

Let us start with a brief formulation of the height 
model for the bipartite N1-N4 dimer model. This is anal- 
ogous to Henley's work on the pure Ni model^ Readers 
interested in deeper detail are referred to Rcfs. 

Crucially, bipartitness allows us to orient each dimer 
so as to point from one sublattice to the other. We can 
then assign to a dimer n units of a fictitious magnetic 
flux, B, and to a link not occupied by a dimer m units 
of flux. The lattice divergence V • B = if we choose 
n + 11m = 0, as in the N1-N4 model, 11 of the 12 links 
emanate from a site are empty. It is convenient to choose 
the overall scale of B so that n = —11m = 11/12. The 



constraint V • B = can be resolved in terms of a scalar 
height function, h, via B — V/i. (In d = 2, curl and 
grad are basically interchangeable). The crucial step in 
the height model ansatz is to posit that, upon coarse- 
graining, the entropy of the height surface gives rise to 
an energy functional containing a leading quadratic term 
of the form 



E = irK d 2 f\Wh\ 



(4) 



Here, the (as yet undetermined) coupling constant K is 
known as the stiffness K. 

Implicit in the derivation of such an energy functional 
of entropic origin is a coarse-graining procedure, which 
maps many microstates of the dimer model onto fewer 
macrostates^ In the process, other terms are generated 
both for the energy and in the identification of operators 
linking the height variable to, say, dimers and monomers. 
These terms are classified by an 'electromagnetic' charge, 
which, together with K, determines the exponent of the 
leading power of its correlationsii&MiSi 

First, let us discuss the N1-N2 model. N2 dimers vi- 
olate bipartiteness and thus show up as defects in the 
height surface; a finite density of such defects destroys 
the algebraic correlations. This was calculated analyti- 
cally for a planar dimer model interpolating between the 
square and the triangular lattices by switching on only 
one of the two diagonal directions of next-nearest neigh- 
bour bonds. There, a dimer correlation length propor- 
tional to the inverse of the small fugacity was also found, 
which result was supplemented by an effective field the- 
ory arguing the same should hold true for monomer cor- 
relations in this system^ Our numerical results, in par- 
ticular the scaling plots, confirm that these results hold 
also for the case of the nonplanar N1-N2 model. 

Next, we consider the N1-N4 model. We first note that 
the monomer correlations have a power-law decay with 
an exponent given by K. The decrease of K from 1/2 
upon adding Ni dimers, implies that the corresponding 
height model grows floppier (the heights fluctuate more 
strongly) . We note in passing that it it is rather straight- 
forward to extract the monomer correlations using the 
directed-loop algorithm, as this only requires keeping 
track of the monomer separations in intermediate states 
sampled. This thus provides a means of determining the 
stiffness constant K with relatively little effort. 

The analytical calculation of the dimer correlations are 
a little more involved than the monomers, as they con- 
tain two terms. One of them, the 'vertex term', decays 
with an exponent of 1/K. The other, the 'dipolar term', 
always decays as 1 /r 2 ; which of these dominates depends 
on K, and for the pure N% model, they happen both to 
decay with the same power law. It is the vertex term 
which generates the logarithmic peak in the D± dimer 
structure factor in Fig. [SJ Upon inclusion of N4 dimers, 
i.e. with increasing K , the power of the decay becomes 
more rapid, and the divergence of this peak thus disap- 
pears when N± dimers are added. 



The survival of the leading 1/r 2 power in the dimer 
correlations is solely due to the dipolar term, which 
arises via the original identification of the dimers and 
V/i through the lattice flux B: If we choose the vec- 
tors ek (k — 1 • • • 12) to point along the direction of a 
dimer from one sublattice to the other, we find a contri- 
bution to the dimer density nk which is proportional to 
efe • B — ik ■ V/i. As long as the height correlations grow 
logarithmically with distance, ((h(r) — h(0)) 2 ) oc lnr, as 
they do for the parameters our simulations probe, the 
resulting contribution to the dimer correlations always 
decay as 1/r 2 , on account of the derivative linking h and 
B. As the power 1/K of the vertex term is only slightly 
larger than 2 for small values of u>4, this leads to the 
very slow approach to the asymptotic 1/r 2 behaviour, as 
evidenced in our simulations. 

Despite their 1/r 2 behaviour, the dipolar correlations 
do not yield a logarithimic peak in the structure factor. 
This is because their dipolar structure makes the lead- 
ing term vanish as the angular part of the Fourier inte- 
gral Dyi(q) — J dr d(fiexp(iq- r) Dyi(r,4>) is performed. 
Instead, what one finds are the characteristic bow-tie 
struct ure o 36 i 37 i 38 displayed in Fig. |S] These bow-ties are 
points of non-analyticity; there is no divergence, but a 
discontinuity there. The bow-ties reflect the transverse 
correlations of a divergence-free field B. These bow-ties 
are thus oriented with one axis along the dimer under 
consideration. 

It has recently been observed 3 ^ that a semi- 
quantitative theory for such correlations can be obtained 
in the framework of a large- N generalisation of the hard- 
core dimer model, via an intermediate mapping to an 
Ising model, the ground-states of which are in correspon- 
dance with the dimer states. This Ising model is then 
solved in a self-consistent mean-field theory in the limit 
of low temperatures. We note that such a generalisation 
does not preserve the charge-assignments for the opera- 
tors generated upon coarse-graining, and does therefore 
miss the vertex term and hence the logarithmic peak for 
the JVi model. However, it does provide the angular de- 
pendence of the dipolar part of the correlations. 

Here, we present the relevant theory for the dimer cor- 
relations both in the pure N\ and the pure N4 models. 
The former is of course also soluble in closed form. In ei- 
ther case, one assigns an Ising spin S to each bond of the 
lattice which encodes the presence (S = +1) or absence 
(S = —1) of a dimer on this bond,3S> The hardcore corre- 
sponds to a non-zero magnetisation in the model, but this 
leaves the shape of its q ^ correlations unchangedi 3 ^ 3 ^ 

In the case of the Ni model, there are two bonds per 
unit cell (the links of the square lattice in the x and y 
directions), denoted by vectors v a , «i,2 = x,y. The rele- 
vant interaction matrix, J, is thence a 2 x 2 matrix. Its 
Fourier transform is given by J a b — cos (q a /2) cos ((ft/2). 
Here, q a — q ■ v a , so thali 3 ^ 




FIG. 7: Fourier transform of the dipolar part (i.e., without the 
dominant logarithmic peaks in Ni) of the dimer correlations 
Di of the Ni model (left) and D4 of the N4 model (right), 
calculated from a large-JV theory. The over-all amplitudes 
are not given accurately by the theory and are therefore not 
indicated here. 



For the N4 model, due to the non-planarity of the lat- 
tice on which the dimer model is defined, this theory 
takes a somewhat more complicated form. One now ob- 
tains a 4 x 4 interaction matrix. Again, let v a denote 
the four different dimer directions: — (2, ±1), u 3)4 = 
(±1,2). The same form of the interaction matrix as above 
then holds: J a b — cos (q a /2) cos (<ft>/2). We thus find the 
following correlations: 



Diiq) ~ 1 



cos 2 (q ■ vi/2) 
ELi cos2 (9- «a/2) 



(6) 



D A (q) 



cos 2 (q y /2) 



cos 2 (q x /2) + cos 2 (q y /2)' 



(5) 



These low-temperature spin (and hence approximate 
dimer) correlations are plotted in Fig. [7| above. We note 
that not only are the bow-ties of the simulation data, 
Fig. [5] reproduced in their correct orientation, but the 
qualitative structure of D4 in reciprocal space is repro- 
duced throughout. 



V. CONCLUSIONS 

The results obtained here demonstrate that the 2D 
square-lattice dimer model becomes deconfined when any 
finite density of non-bipartite (next-nearest-neighbor) 
dimers are introduced. The system remains alge- 
braically confined in the presence of longer bipartite 
(fourth-nearest-neighbor) dimers, but the corresponding 
monomer exponent drifts. In contrast, as the stiffness of 
the corresponding height surface decreases, the leading 
dimer-dimer correlation exponent remains unchanged at 
1/r 2 , due to the persistance of the dipolar term due to 
the gauge structure captured by the height model. Our 
results thus fit in with the known beliefs on dimer mod- 
els and can be seen as a detailed check of their validity. 
We have shown that the directed-loop simulation algo- 
rithm can produce very accurate results for the drifting 
monomer exponent. In particular, the accuracy is high 
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enough to suggest the conjecture that this exponent for 
the pure N4 model takes the rational value of 1/9. 

Our findings are relevant to quantum dimer models as 
well. Various resonance terms can be introduced, and 
corresponding RK points can then be demonstrated in 
the same way as has been done for other models ^! 14 - 4 ? 
Hence, it is clear that an RVB state can in principle be 
realized on the square lattice once next-nearest-neighbor 
dimers are allowed. In the context of the dimer models 
for the Mott instulators relevant to the high-T c cuprates, 
two points however need to be borne in mind. Firstly, the 
Rokhsar-Kivelson quantum dimer models allows a simple 
classical calculation at the RK point only if it is possible 
to choose the sign of the off-diagonal matrix elements to 
be negative, which property needs to be established on a 
case-by-case basis for non-bipartite models (it does hold 
for the bipartite case)4i In addition, the description of 
the singlet space using a dimer basis becomes overcom- 
plete as the coordination of the dimer lattice increases. 
At any rate, our results again underscore the special im- 



portance played by non-bipartitness in the realisation of 
RVB liquids in d = 2 + 1. 
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